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A strongly well-posed initial boundary value problem based upon constraint-preserving bound- 
ary conditions of the Sommerfeld type has been established for the harmonic formulation of the 
vacuum Einstein's equations. These Sommerfeld conditions have been previously presented in a 4- 
dimensional geometric form. Here we recast the associated boundary data as 3-dimensional tensor 
fields intrinsic to the boundary. This provides a geometric presentation of the boundary data analo- 
gous to the 3-dimensional presentation of Cauchy data in terms of 3-metric and extrinsic curvature. 
In particular, diffeomorphisms of the boundary data lead to vacuum spacetimes with isometric ge- 
ometries. The proof of well-posedness is valid for the harmonic formulation and its generalizations. 
The Sommerfeld conditions can be directly applied to existing harmonic codes which have been used 
^ ' in simulating binary black holes, thus ensuring boundary stability of the underlying analytic system. 

O ' The geometric form of the boundary conditions also allows them to be formally applied to any metric 

, formulation of Einstein's equations, although well-posedness of the boundary problem is no longer 

ensured. We discuss to what extent such a formal application might be implemented in a constraint 
(N ■ preserving manner to 3-1-1 formulations, such as the Baumgarte-Shapiro-Shibata-Nakamura system 

which has been highly successful in binary black hole simulation. 
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I. INTRODUCTION 



I Previous work has established the strong well-posedness of the initial-boundary value problem (IBVP) for Einstein's 
^ • equations expressed as a hyperbolic system in harmonic coordinates. The result was first obtained using pseudo- 

' differential theory, i.e. using a Fourier-Laplace expansion, which established well-posedness in the generalized sense 
00 : The result was subsequently obtained using standard energy estimates 0, i]. This places the IBVP on the same 
0^ analytic footing as the Cauchy problem, whose well-posedness was also established using harmonic coordinates in 
' I , the classic work of Choquet-Bruhat Q. The geometric formulation of the boundary conditions and boundary data 
0^ ' for the IBVP is more complicated than for the Cauchy problem. Recently, this boundary data was presented in a 
' 4-dimensional geometric form Q. Here we recast the boundary data as 3-dimensional tensor fields intrinsic to the 
boundary, analogous to the presentation of Cauchy data in terms of the 3-metric and extrinsic curvature of the initial 
Cauchy hypersurface. The spacetime metric which solves the harmonic IBVP with this data is uniquely determined 
up to a diffeomorphism. 

In the Cauchy problem, initial data on a spacelike hypersurface Sq determine a solution in the future domain of 
dependence T>^{So) (which consists of those points whose past directed characteristics all intersect Sq). In the IBVP, 
■ data on a timelike boundary T which meets Sq in a surface Bq are used to further extend the solution to the domain of 
dependence V^(Sq UT). In practical applications Bq is topologically a sphere surrounding some system of interest but 
here we only consider the local problem in a neighborhood of a point of intersection between the Cauchy hypersurface 
and the boundary. For hyperbolic systems, the global solution in the spacetime manifold A4 can be obtained by 
patching together local solutions. The setting for this local problem is depicted in Fig. [D 

There are no natural boundaries for the gravitational field analogous to the conducting boundaries that play a major 
role in electromagnetism. Consequently, the IBVP for Einstein's equations only received widespread attention after 
its importance to the artificial outer boundaries used in numerical relativity was pointed out [Gj] . The first well-posed 
IBVP was achieved for a tetrad formulation of Einstein's theory in terms of a first differential order system which 
included the tetrad, the connection and the curvature tensor as evolution fields f^. A strongly well-posed IBVP was 
later established for the harmonic formulation of Einstein's equations as a system of second order quasilinear wave 
equations for the metric [l|, Q ■ Strong well-posedness guarantees the existence of a unique solution which depends 
continuously on both the Cauchy data and the boundary data. The results were further generalized in 3] to apply to 
general quasilinear symmetric hyperbolic systems whose boundary conditions have a certain hierarchical form. 

The initial data for the Cauchy problem can be formulated in terms of two symmetric tensor fields hat (a Euclidean 
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FIG. 1; Data on the 3-manifolds So and T, which intersect in the 2-surface Bo, locaUy determine a solution in the spacetime 
manifold M 



S-metric) and kat on a 3-manifold Sq subject to the (Hamiltonian and momentum) constraints 
and 

= Vaikb-Sfk'J, (1.2) 

where Va is the covariant derivative and ^^^i? is the curvature scalar associated with hat- As characterized in [9*], via 
an embedding in a 4-dimensional manifold A4, the data hab and kab on the disembodied 3-manifold Sq determine a 
Lorentzian metric satisfying Einstein's equations which is unique up to diffeomorphism and whose restriction to the 
embedding Sq oi Sq gives rise to its intrinsic 3-metric hab and extrinsic curvature kab- 

Here we present an analogous result for the IBVP. The difhculties underlying the IBVP, which have most recently 
been discussed in Q and Q , make the formulation of disembodied boundary data more difficult than for the Cauchy 
problem. There are three main complications. 

• The first complication stems from a well-known property of the IBVP for the flat-space scalar wave equation 

{df - V^)* = 

in the region a; < 0. Although the initial Cauchy data consist of ^\t=o and dt^\t=o, only half as much boundary 
data can be freely prescribed at a; = 0, e.g the Dirichlet data dt^\x=o, or the Neumann data dx^\x=o or the 
Sommerfeld data {dt+dx)^\x=o (based upon the derivative in the outgoing characteristic direction). For a given 
physical problem, this implies that the boundary data cannot be prescribed before the boundary condition is 
specified, i.e. the correct boundary data depend upon the boundary condition, unlike the situation for the 
Cauchy problem. The analogue in the gravitational case is the inability to prescribe both the metric and its 
normal derivative on a timelike boundary, which implies that you cannot freely prescribe both the intrinsic 
metric of the boundary and its extrinsic curvature. This leads to a further complication regarding constraint 
enforcement on the boundary, i.e. the Hamiltonian and momentum constraints (|l.ip and (|1.2p cannot be enforced 
directly because they couple the metric and its normal derivative. Here we restrict our attention to Sommerfeld 
boundary data which we prescribe in a constraint free manner. 

• A Sommerfeld boundary condition for a metric component supplies the value of the derivative K'^dcQab in an 
outgoing null direction K'^. Such boundary conditions are very beneficial for numerical work since they allow 
discretization error to propagate across the boundary (whereas Dirichlet and Neumann boundary conditions 
reflect the error and trap it in the numerical grid). However, the boundary does not pick out a unique outgoing 
null direction at a given point (but, instead, essentially a half null cone). This complicates the geometric for- 
mulation of a Sommerfeld boundary condition. In addition, constraint preservation does not allow speciflcation 
of Sommerfeld data for all components of the metric, as will be seen in formulating the Sommerfeld conditions 

- m^- 
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• The third major comphcation arises from gauge freedom. In the evolution of the Cauchy data it is necessary to 
introduce a fohation of the spacetime by Cauchy hypersurfaces St , with unit timchkc normal Ua ■ The evolution 
of the spacetime metric 

gab = -naUf, + hab (1-3) 

is carried out along the flow of an evolution vector field t"' related to the normal by the lapse a and shift /?" 
according to 

t" = an'' + , (3''na = 0. (1.4) 

The choice of foliation is part of the gauge freedom in the resulting solution but does not enter into the 
specification of the initial data. In the IBVP, the foliation is unavoidably coupled with the formulation of the 
boundary condition. Thus some gauge information must be incorporated in the formulation of the boundary 
condition and boundary data. 

In order to resolve these complications, we include the specification of a foliation Bt of the boundary T as part 
of the boundary data. This supplies the gauge information which determines a unique outgoing null direction for 
a Sommerfeld condition. In Section |TT1 we formulate the 3-dimensional prescription of the boundary data and the 
statement of our main result, a Theorem establishing the existence of a solution satisfying a version of the geometric 
uniqueness property proposed in [8|]. The 3-dimensional description of the data in terms of scalar, vector and tensor 
fields intrinsic to the boundary is quite abstract at this stage. Their physical significance only becomes clear in 
Section [mi where the proof of the Theorem is given. The proof is based upon the well-posedness of the 4-dimensional 
geometrical formulation of the harmonic IBVP given in [5|. Constraint preservation is established by incorporating 
the harmonic conditions into the boundary conditions. 

The motivation for this work stems from the need for an improved understanding and implementation of boundary 
conditions in the computational codes being used to simulate binary black holes. We discuss the applicability of our 
results to numerical relativity in Sec. IIVI 

Much of the presentation in the paper is coordinate independent and we use Latin letters (a, 6, c, ...) as abstract 
indices [lol | to denote the types of vector and tensor fields and to indicate their manipulations. This notation serves 
to describe either 4-dimensional tensor fields on the spacetime M or 3-dimensional tensor fields intrinsic to Sq or T. 
When spacetime coordinates = {t, x*) are introduced, we use Greek letters (/i, ly, p, ...) to describe the corresponding 
4-dimensional tensor components and Latin letters (i, j, k, ...) to denote the spatial components. 

II. DISEMBODIED DATA FOR THE IBVP 

We state our main result concerning Sommerfeld boundary data for Einstein's equations. 

Geometric Uniqueness Theorem: Consider the 3-manifolds T and Sq meeting in an edge Bq. On Sq prescribe 
the smooth, symmetric tensor fields hab and kab, subject to the Hamiltonian and momentum constraints and the 
condition that hab be a Riemannian metric . On Bq prescribe the smooth scalar field O. On T prescribe a smooth 
foliation Bt parametrized by a scalar function i, where (for convenience) t = on Bq. In addition, on T, prescribe 
the scalar field q, the vector field q"" and the rank-2 symmetric tensor field a"^, which are all smooth and vanish on 
Bq- Here the rank-2 property of a"'' is defined with respect to the foliation Bt by the requirement 

a°'''dai=0. (2.1) 

Then, after embedding 5o U T as the boundary iSq U T of a 4-manifold M, as depicted in Fig. HI this data provides 
Sommerfeld boundary data for a vacuum spacetime, in a region including a neighborhood of the embedded edge Bq, 
which is unique up to diffeomorphism. 

The geometrical interpretation of the data involves the metric of the embedded spacetime, whose existence is the 
content of the Theorem. Before proceeding to the proof in the next section, it is useful to supply some intuitive 
meaning to the data. As in the Cauchy problem, hab and kab are identified with the 3-metric and the extrinsic 
curvature of the embedding So of So- The scalar field O determines the hyperbolic angle describing the initial velocity 
of the embedded boundary T relative to the inertial frame picked out by Sq. The fields q and q°' supply information 
concerning the subsequent dynamics of the boundary and its foliation. Together q and determine the components 
of a 4- vector describing the curvature of the outgoing null geodesies normal in A4 to the embedding Bt of Bt- The 
field a"'^ determines the optical shear a of the outgoing null hypersurface through Bt- 
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Note that the data contains no metric information about the embedded boundary T, not even that it is a timehke 
3-manifold. Such structure only emerges from the construction of a solution for the spacetime metric. Thus fields 
derived algebraically from the metric, such as the unit normal to the boundary, are to be considered as subsidiary 
unknowns. 

The requirement that the boundary data vanish at Bq stems from its interpretation relative to the initial Cauchy 
data. This requirement ensures the continuity of the resulting spacetime metric and its first derivatives. The full 
compatibility conditions necessary for a C°° spacetime metric involve enforcing the Einstein equations and their 
derivatives on Bq. It is not clear how to implement C°° compatibility conditions in terms of disembodied data. 



III. THE EXISTENCE OF A GEOMETRICALLY UNIQUE SOLUTION 

Our task is to show that harmonic evolution of the data prescribed in the Geometric Uniqueness Theorem of Sec. |TT] 
locally determine a vacuum spacetime which is unique up to diffeomorphism. We embed Sq, T and Bq in a 4- manifold 
A4 with boundary consisting of the corresponding pieces Sq, T and Bo. The unknown is a spacetime metric gab on M 
which satisfies Einstein's equations and is determined up to a diffeomorphism by the data specified in the Theorem. 
We follow the construction given in Q. 

The first step is to give initial data for gab- Let hab and kab be the fields induced by hab and kab on iSq; Q be the 
field induced by 6 on i3o; and g, q"' and (T°^ be the fields induced by q, q"" and ct"'' on T. Let Bt be the fohation of 
T corresponding to the embedding of Bt, where t is the parametrization induced by i. On Sq prescribe a transverse 
field and construct the Lorentzian metric g"''' — —n°'n'' + /i"^, so that is the future directed unit normal to Sq. 
Require that n° satisfy n'^Naleo — sinhO, where Na is the outward unit normal to T. 

We introduce the boundary decomposition of the metric 

gab^ NaNh~TaTb + Qab, (3.1) 

where Ta is the future directed unit normal in T to Bt- (Here the boundary fields Na, Ta and Qab are unknowns 
subsidiary to gab-) This leads to an orthonormal tetrad (T", A^", Q", Q") on T, where is a complex null vector 
tangent to Bt with normalization 

Qab = QiaQb), Q"Qa = 2, Q'^Qa^O- (3.2) 

(The tetrad is unique up to the spin freedom — > e^^Q" which docs not enter our construction in any essential way.) 
Uniquely associated with this tetrad (independent of the choice of Q'^) are the outgoing and ingoing null vector fields 
j^a _ rpa _^ jya ^^^^ L"- = T"" — N"" , respectively, which lie in the null directions normal to Bt- They form a null tetrad 
{K°- , L°',Q'^ ,Q°-) with metric decomposition 

gab = -K(aLb) + Q(aQb)- (3.3) 

Next we tie down the gauge freedom by introducing an evolution field on We require that t"" be tangent to 
T and that it generate the foliation Bt according to 

Ctt = I, (3.4) 

where Ct is the Lie derivative with respect to t"' - We extend to AA such that it generates a foliation St, with 
parametrization satisfying (|3.4|) and unit normal n°. The gauge freedom is then pinned down by introducing spatial 
coordinates a;* on So and extending them to M according to 

Ctx' = 0. (3.5) 

The scalars a;^ = {t, x^) serve as coordinates for M which are adapted to the evolution. Note that <° and the adapted 
coordinates x^^ are explicitly constructed fields on with no metric properties. In the relationship (|1.4p between 
and the unit normal n° to St, it is n° which contains metric information and is a subsidiary unknown. 

As part of the initial data, we prescribe Ctn"" on Sq- Along with the initial choice of n°, hab and kab, this determines 
the initial data 5afc|t=o and Ltgab\t=o- We use the evolution field to provide a background metric (J^f, on AA which 
is uniquely and geometrically determined by the Lie transport of the initial data according to 

gab\t=o = gab\t=o , Ctgab\t=a = Ctgab\t=a , ^tCtgab = 0- (3-6) 
In the coordinates x'^ = {t,x^) adapted to the evolution, this reduces to 



(3.7) 
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Ct^ = ( Va/bc + ^ahc - VJab ) . (3.10) 



The connection Vi a-nd curvature tensor R'^cab associated with the background g^f^ have the same transformation 
properties as the corresponding quantities Vq and R'^cab associated with gab- In particular, the difference Va — 
defines a tensor field Cf^ according to 

(Va - = Gtbv\ (3.8) 

for any vector field . In terms of the (nonlinear) perturbation 

fab = gab ~ 9ab (3-9) 

of the metric from the background, we have 

1 

T 

We take fab to be the evolution variable for solving Einstein's equations. Since gab is explicitly known, a solution 
for fab is equivalent to a solution for gab- By construction, the initial data for fab is homogeneous, i.e. 

/afc|t=0 = Ctfablt^O = 0. (3.11) 

The boundary data on T consist of the vector field 

q'' = qN'' + q'' (3.12) 

and the tensor field 

a"" = - ^Q'^'Qcd^^''. (3.13) 

Here, and elsewhere, the physical metric gab is used to raise and lower indices and to normalize the tetrad vectors. 

Now, just as the Cauchy data hab and kab must be identified as the intrinsic metric and extrinsic curvature of Sq 
in order to construct a solution of the initial value problem, we give a geometric identification of the boundary data. 
We identify this data as the geodesic curvature and shear, relative to their background metric values, of the outgoing 
null vector K"" on T according to the formulae 

= K\\7b~^)K\ (3.14) 

= ^(g'^^Q'"'-^Q'^V)(Vc-V,)ifd. (3.15) 
The last equation can be re-expressed in the spin-weight-2 form 

a Q'^gVafc - ^Q°Q'(Va - Va)Kb- (3.16) 

The use of the shear in posing geometrical boundary conditions for the harmonic formulation was suggested earlier 
in 

Using p.Sp and p.lOp . we recast p.l4p and (|3.16p as Sommerfeld boundary conditions which determine the com- 
ponents of the outgoing null derivatives K'^Vafbc according to 

Ix'-K-K'^Vafbc - g"i^a, (3.17) 

{Q^'K-K'^ - \K''K-Q'^)Vafbc - g"Qa, (3.18) 



2 

[L^'K^K'^ ^]^K'-K-L'^)Vafbc = q'^La, (3.19) 



(^Q'Q'if" - Q'^'Q")Va/bc = 2f7. (3.20) 

In addition to these six Sommerfeld conditions, we impose the four additional boundary conditions that :— 
gabfjd^ = on T, i.e. the harmonic constraints. In terms of the null tetrad decomposition, they take the Sommerfeld 
form 

- 2C''Ka = [Q^Q^K" + K''K''L'' - K'^Q'^Q" - K^Q^Q") Vafbc ^ 0, (3.21) 
-2C"ga = (i'-g^i^" + K^'Q^L'^ - K'^L^Q'^ + Q^'Q^Q'^) Vafbc - 0, (3.22) 
-2CUa = iL^'L^K" + Q^'Q^L" - Q^'L^Q" - Q'^L^Q'^) Vafbc = 0- (3.23) 
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Together, (|3.17p - (I3.23P provide Sommerfeld boundary conditions for the components of K°-'SIafbc in the sequential 
order (KK), {QK), {LK), (QQ), {QQ), {LQ), {LL) in terms of the boundary data and the derivatives of preceding 
components in the sequence. Such a hierarchy of Sommerfeld boundary conditions satisfy the requirements of Theorem 
1 of [3] which establishes a strongly well-posed IBVP for a quasilinear hyperbolic system. 

In order to apply this theorem, we reduce Einstein's equations to a quasilinear wave system by modifying the 
Einstein tensor by the harmonic constraints 

C = g'^'C^b = Ig'^'g'' {^ahc + ^ahc - kja^ (3.24) 

according to 

^ab_^ab__y(a^b)^ l^afcy^^rf^ (3.25) 

We eliminate the coordinate freedom, up to our choice of evolution field t° on M. and coordinates on iSq, by working 
in the coordinates adapted to the evolution. The coordinate components of p.24|) then take the form 

C = ig'^^s"" (i/.. + - = s'^^^U - h^)- (3.26) 

The requirement that V :— = is the standard harmonic coordinate condition that Dga;^ = 0. In the present 

case, setting C = implies — g'^'^T'^^, which is an example of harmonic coordinates with a forcing term, as 
discussed in T^]. (In numerical relativity, these have been called generalized harmonic coordinates [3.) For more 
general forcing terms the harmonic constraint takes the form 

V^HPix^g), (3.27) 

where restriction of the forcing term to depend only upon x^'and g^^ ensures that the system remains well-posed. 
In the adapted coordinates, the reduced Einstein equations i^^jy = form the desired quasilinear wave system for 

S^'ii/M- = ^gxrg'^C^pC^^ + 4C^(^g,)^C%,5"" - 25''"i\.(^3.)A. (3.28) 

With the hierarchy of Sommerfeld boundary conditions (|3.17p - p.23p , Theorem 1 of [3] now applies and ensures that 
(I3.28P has a well-posed IBVP and, in particular, determines a unique solution f^^. In addition, the resulting metric 
(7^^ must satisfy the harmonic constraints C = because they are built into the initial data and boundary conditions. 
(See Sec. llVl for details concerning constraint preservation.) Therefore g^i^ solves Einstein's equations. 

The solution has been obtained in coordinates which are harmonic with respect to the background g^,^, i.e. 5^''(r^^ — 

r^y) = 0. The resulting spacetime metric gab and background metric g^f, determined by the data are in a gauge which 
depends upon the choice of evolution field t"" . Under a diffeomorphism \1/ of which reduces to the identity map 
on So and T, we have {f^ , gab, gab) ~^ ^*5ab> ^*5a&)- Thus 'i'*gab is a diffeomorphic solution of Einstein's 

equations which is harmonic with respect to the background gat,- 

A more important question, which concerns the issue of geometric uniqueness raised in [8| , is the behavior of the 
spacetime under a diffeomorphism 5' of the disembodied boundary T. It is known for the Cauchy problem that 
data {^*hab, ^*kab) on So lead to a spacetime which is isometric to the spacetime with Cauchy data {hab, kab)- The 
analogous result holds for the boundary data. Under a diffeomorphism ^ of T the boundary data maps according 
to (f, g, CT°^) ('l*^,^*g,^,g°^^',a•°''). After the embedding in M, this induces data (5'*i, *,(t''^) 
on T. Now consider any smooth extension of 4* to a diffeomorphism of A4. Suppose the original data leads by 
the above construction to the spacetime metric gab with the evolution field background gah and relative geodesic 
curvature and relative shear {q°',a°'^) of the outgoing null vector if normal to the foliation of T. Then, by the same 
construction, the mapped data lead to the metric ^*gab with the evolution field vE'*i°, background ^*gab and relative 
geodesic curvature and relative shear (^^g", ^*ct°^) of the outgoing null vector \P,i4r" normal to the ^'*i-foliation of 
T. In this way, diffeomorphisms of S and T, which map their intersection Bo into itself, generate an isometry class 
of vacuum spacetimes. 

We have thus shown that the resulting vacuum spacetimes are diffeomorphic if the disembodied boundary data and 
initial Cauchy data are diffeomorphic. This comprises a version of geometric uniqueness, which has been pointed out 
as a missing ingredient in prior formulations of the IBVP However, the result is not as strong as for the pure 
Cauchy problem for which the converse is also true: the resulting vacuum spacetimes are diffeomorphic only if the 
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initial data are difFeomorphic. The converse is more complicated for the IBVP and does not hold in the disembodied 
setting because the data a""^ are superfluous. It is only the trace-free part of the embedded data ct'^'' which enters the 
shear a°'^ . But there is no way to pick out the trace free part without knowledge of the metric, which is an unknown 
at the stage of specifying data. (For linear perturbations, a background metric could be used but this construction 
does not extend to the nonlinear case.) This is perhaps an unavoidable feature of disembodied Sommerfeld data for 
the IBVP. 

As a result, disembodied boundary data which are not related by a diffeomorphism can lead to isometric spacetimes. 
However, it follows by direct geometrical construction that an isometry class of spacetimes does uniquely determine 
the embedded boundary data and ct"'' up to a diffeomorphism. It is only upon the transition from ct"'' to (t°^ 
that the ambiguity of adding a trace enters. In this sense, a stronger version of geometric uniqueness applies to the 
embedded data. 

The Sommerfeld boundary conditions were based upon the background metric obtained by the Lie transport (j3.6p 
of the Cauchy data along the streamlines of i". Modification of this transport law would lead to a different background 
and the spacetime generated by the same disembodied data would in general not be diffeomorphic. This is similar 
to the scalar wave problem, for which the same boundary data would lead to different solutions if, say, a Dirichlet 
boundary condition were used instead of a Sommerfeld condition. In the present case, a new background metric 
implies a different form of the Sommerfeld condition and, for fixed boundary data, the resulting spacetime will (in 
general) not be isometric to the original spacetime. The effect of the boundary data is changed by the difference 
between the two background connections. 

IV. NUMERICAL APPLICATION 

The motivation for this work has been the treatment of the outer boundary in the numerical simulation of the 
inspiral and merger of binary black holes. The boundary conditions (|3.17p - (|3.23p . which lead to a strongly well- 
posed IBVP, can be applie d directly to any of the harmonic evolution codes which have been used to simulate this 
binary problem [T^. [Tsl. ITgI. [TtI. [Tsj . At present, none of these harmonic codes incorporate boundary conditions that 
ensure strong well-posedness. The closest example is the pseudo-spectral harmonic code described in [l^ [l3| which 
incorporates a second differential order boundary condition which freezes the Weyl component and was shown to 
be well-posed in the generalized sense in the high frequency limit JJj]. 

An important attribute of strong well-posedness is the estimate of the boundary values of the solution and its 
derivatives which are provided by an energy conservation law obeyed by the principle part of the equations. This 
boundary stability extends to the semi-discrete system of ordinary differential equations in time which are obtained by 
replacing spatial derivatives by finite differences obeying summation by parts (the discrete counterpart of integration 
by parts), so that energy conservation caries over to the semi-discrete problem. This stability then extends to the 
numerical evolution algorithm obtained by applying an appropriate time integrator, such as Runge-Kutta. 

The advantage of boundary stability in numerical applications is that the smoothness of the solution is unaffected 
by reflections off the boundary. This avoids the long timescale instabilities which might otherwise arise from multiple 
reflections. For the harmonic Einstein problem, each component of the metric obeys a quasilinear wave equation so 
that it is straightforward to develop a summation by parts algorithm based upon the standard energy expression for 
a scalar wave in a curved spacetime ^19] . A code incorporating such an algorithm was applied, using a version of 
the well-posed Sommerfeld boundary conditions presented here, to the test problem of a highly nonlinear gauge wave 
propagating inside a cubic boundary . Although the proof of strong well-posedness given in [3| was based upon a 
scalar wave energy differing by a small boost from the standard energy expression, the successful results for this test 
problem confirm the robustness of the underlying approach. 

Strong well-posedness extends to more general quasi-harmonic formulations for which the reduced Einstein equations 
have the form 

^M"^ Qi^i' „ v'^C") -I- iff'^^VpC -f A^'^C" = 0, (4.1) 

where C = g'^T^^ — HP{x,g) are the generalized harmonic constraints (|3.27p and the coefficients A^'^ have the de- 
pendence A^'^(x, g, dg). (This includes constraint modified versions of the harmonic system.) Constraint preservation 
follows from the Bianchi identity V ^G^'^ = which implies a homogeneous wave equation for C, 

VPVpC + Ri'pCP ~2Vp{A^PC'') ^0. (4.2) 

If the boundary conditions enforce Cp\t = and the initial data enforces C^Isq = dtCP\so — then the unique solution 
of (|4.2p is C = 0. As a result, the Sommerfeld boundary conditions in the geometrical form (|3.17p - p.20p . along 
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with l|3.2ip - (|3.23p which enforce Clr = 0, lead to a well-posed harmonic IBVP in which the harmonic constraints 
C = are satisfied everywhere. In turn, (|4.ip implies that the Hamiltonian and momentum constraints G^'^rii, — 
are also satisfied. 

Beyond the geometrical aspects of the harmonic IBVP, there are practical concerns that arise in astrophysical 
applications. The linear time dependence of the background metric p.7p could lead to deleterious long time scale 
effects in numerical simulations. For that reason, it is preferable to modify the prescription (13. 6p for the background 
metric so that (13.71) changes to 

9^,u=9^^u\t=o+te-^'{^tg^J.^t=o) (4.3) 

and the time dependence damps on a time scale determined by A. Alternatively, the purpose of tying the background 
metric to the initial Cauchy data was to make clear that it did not affect the geometric uniqueness of the solution. 
Otherwise, the Minkowski metric in the coordinates adapted to the evolution could have been chosen as the background 
metric. This would leave intact the well-posedness of the IBVP and might be the most expedient approach in some 
numerical applications, although (j4.3p has the advantage of suppressing nonlinear effects in the early stage of an 
evolution. 

Another practical concern is that the proper boundary data and a are not normally known in an astrophysical 
application. The practice in simulating an isolated system is to assume homogeneous data on the artificial outer 
boundary, i.e. = cr = 0. Such homogeneous data in general leads to some spurious back reflection from the 
boundary. However, as discussed in fsl for the harmonic IBVP, similar boundary conditions on a round spherical 
outer boundary of large surface area radius R lead to reflected waves whose amplitude falls off asymptotically as 1 /B? 
for both quadrupole gauge waves and quadrupole gravitational waves. (The calculation assumes that the linearized 
approximation is applicable. Modifications of the boundary conditions by lower differential order terms involving 
factors of R lead to a faster falloff for the reflected quadrupole gravitational waves. 

The boundary conditions (|3.17p - (|3.23p are slightly different than those considered in 3] and, in particular, the 
homogeneous boundary condition a — Q leads to reflected waves with only a falloff, as previously noted in [ll[ . 

However, improved performance can be obtained by taking advantage of the sensitivity of the boundary conditions 
to the location of the indices. If p.l6p and (|3.20p are replaced by 

then agreement with the corresponding boundary condition (94) of 0] is obtained. In that case, for the boundary 
conditions (|3.17p - (|3.19p and (|4.4p , along with the harmonic constraints (|3.2ip - (|3.23p , the amplitudes of the reflected 
quadrupole gauge waves and quadrupole gravitational waves again fall off asymptotically as Thus the application 

of these Sommerfeld boundary conditions with homogeneous data results in small back reflection from the boundary 
of an isolated system. 

While the Sommerfeld boundary conditions considered here were developed for the harmonic IBVP, their geometric 
nature allows them to be formally applied to any metric version of the reduced Einstein equations, in particular 
the alternative "3 -f 1" formulations upon which much numerical work has been based. However, in that case, well- 
posedness and constraint preservation do not necessarily follow. One approach to dealing with these issues would be 
to re-express the 3 + 1 formulations in terms of the covariant 4-dimensional Z4 theory [28l . [29| , in which the reduced 
evolution equations take the form (|4.ip but with the harmonic constraints replaced by a more general vector fleld 
Z'^, which can be used to change the principle part of the resulting evolution system. The ZA formalism has been 
shown to encompass the standard 3 -f 1 formulations, including the Arnowitt-Deser-Misner (ADM) [s^l , the Kidder- 
Scheel-Teukolsky [3l|, the Nagy-Ortiz-Rcula [s?] and the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) [2l|, ^ 
formulations. It is possible that the close analogue between the ZA and harmonic formulations might be used to 
shed light on the analytic properties of 3 -I- 1 systems. However, such an investigation of the well-posedenss of the 
3-1-1 IBVP is stymied by the fact that none of the systems used in numerical applications have been shown to be 
symmetric hyperbolic. Some 3 -f 1 systems have been show to be strongly hyperbolic, which ensures a well-posed 
Cauchy problem, but symmetric hyperbolicity is required to apply the standard theorems for a well-posed IBVP. For 
that reason, the discussion here is restricted to how the Sommerfeld boundary conditions developed here might be 
applied in a manner consistent with what is presently known about 3-1-1 formulations. 

The boundary conditions ((3l7)) - ((3?20)) can be translated in a straightforward way into a "3-1-1" decomposition 
= (t, x^) = (i, X, y.z), in which the metric takes the form (|1.3p where, in terms of the lapse and shift (|1.4p . 

gtt = -a^ + hijfS'^P^ , gt.^^hijf3\ g^j^hij. (4.5) 

Although there is a different decomposition (j3.ip intrinsic to the boundary, the outgoing null direction used in 
formulating a Sommerfeld condition is intrinsic to both decompositions. Let iV be the unit normal to the boundary 
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foliation Bt, which lies in the t = const Cauchy hypersurface, so that the metric (|1.3p has the further decomposition 

gf^v = -ni_,n^ + N^Ni, + Q^^, (4.6) 

where again Q^jy = Q(f_iQv) is the 2-metric intrinsic to the foliation Bt of the boundary. Then the outgoing null vector 
K^^ normal to the foliation has components K^^ = + = e^^{n^ + N^) where Q is the hyperbolic angle arising 
from the velocity of the boundary relative to the Cauchy hypersurfaces. In terms of the lapse and normal component 
of the shift, 

tanhe = --^*7Vi, (4.7) 
a 

where P^Ni < for a boundary which is moving inward with respect to the Cauchy hypersurfaces. The sign of P^Ni 
determines whether the advective derivative n^9^ is outward < 0) or inward {(3^Ni > 0) at the boundary and, 

as a result, can affect the number of required boundary conditions. 
The Sommerfeld derivative takes the simple 3+1 form 

K'^d^ = e-®(7i^a^ + iv'^a^), = N'd,, (4.8) 

where rt^9^ is the part containing the time derivative dt- Thus the Sommerfeld boundary conditions (|3.17[) - (|3.20[) 
can be used to supply boundary values for the time derivatives of 6 metric components, or equivalently to supply 
boundary values for 5 components of the extrinsic curvature fc^^ of the Cauchy foliation and for the time derivative 
of one metric component. In particular, (|3.17p - (|3.18p and (|3.20p supply boundary values for N^N^k^^, Q^TV'^fc^iy 
and Q^Q'^kfj^^ and (|3.19p supplies the boundary value for the time derivative of the normal component of the shift 
NidtP^, i.e. (/i^^)~^/^9t/3^ for a boundary aligned with the z-coordinate. 

The remaining Sommerfeld boundary conditions (|3.2ip - (|3.23p . which enforce the harmonic constraints, require 
modification depending upon the particular 3 + 1 formulation and gauge conditions. Only 6 components of Einstein's 
equations are used in a 3 + 1 evolution, with the gauge conditions determining the evolution of the lapse and shift. 
As a first example, consider the ADM formulation in which the 6 Einstein equations 

hP^KRp, = (4.9) 

are evolved. The evolution of the Hamiltonian and momentum constraints H := G^i^n^ny and := h^'^rCG^^ 
is governed by the contracted Bianchi identity V^GJ;^ = 0, which gives rise to the symmetric hyperbolic constraint 
propagation system 

n'^d.yH-djP^ = B'^G^jTi" (4.10) 
rfd^P' -K'^djH = B'^^G^^n", (4.11) 

where the coefficients B^ and B^^^ arise from Christoffel symbols and do not enter the principle part. An analysis 
of this system shows that only one boundary condition is allowed provided (3^Ni < 0, i.e provided the boundary is 
moving inward relative to the Cauchy hypersurfaces. The theory of symmetric hyperbolic systems then guarantees 
that all the constraints will be preserved if 

H + P'Ni = G^^n'^K" = (4.12) 

is satisfied at the boundary. (Additional boundary conditions are necessary for constraint preservation if > 0.) 
By combining the evolution system (|4.9p with (|4.12p . this boundary condition is equivalent to 



i.e. the Raychaudhuri equation (cf. [Tc 



G^.K^K'' = 0, (4.13) 



K^'df,e + -e-" + a(f = 0, (4.14) 



where 6 — Q^^'^V^K^, is the expansion of the outgoing null rays tangent to K^. Thus, for the ADM system, constraint 
preservation can be enforced by the Sommerfeld boundary condition (|4.14p for 6, which supplies the boundary values 
for the remaining component Q^'^k^^ of the extrinsic curvature. Unfortunately, although the subsidiary constraint 
system is symmetric hyperbolic, the ADM evolution system is only weakly hyperbolic and consequently leads to 
unstable evolution. 
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In terms of astrophysical applications, the most important 3 + 1 formulation is the BSSN system, which has been 
used by the majority of groups [H, [l^, [H, [H, [23| carrying out binary black hole and neutron star simulations. The 
development of the BSSN formulation has proceeded through an interplay between educated guesses and feedback 
from code performance. Only in hindsight has its success spurred mathematical analysis, which has shown that 
certain versions are strongly hyperbolic and thus have a well-posed Cauchy problem 32^ , ,33 , ^4] . Although significant 
progress has been made in establishing some of the necessary conditions for wcU-posedness and constraint preservation 
of the IB VP m, [1^, [13, nil , there is still no satisfactory mathematical theory on which to base numerical work. In 
current numerical practice, the boundary conditions for BSSN evolution systems are applied in a naive, homogeneous 
Sommerfeld form to each evolution variable (cf. [27|). 

The geometric nature of the Sommerfeld boundary conditions p.l7p - (j3.20p and their role in a well-posed harmonic 
IB VP suggest that they might lead to improved performance over the present boundary treatment of the BSSN system. 
However, there are two complications. The first involves the sign of the normal component of the shift (i'^Ni at the 
boundary, which also entered the above discussion of the ADM constraint system. A recent analysis [11] of the BSSN 
evolution system shows that the number of incoming fields at the boundary, and therefore the number of required 
boundary conditions, depends upon the sign of (3^Ni. A practical scheme for dealing with this would require the 
Dirichlet boundary condition (3^Ni = (or some similar Dirichlet condition to control the sign) rather than the 
Sommerfeld condition (|3.19p for the normal component of the shift. 

The other complication for the BSSN system involves constraint preservation. The BSSN evolution system enforces 
the 6 Einstein equations 

hP^KRp, - ^h^,H = 0, (4.15) 

for which the constraint system implied by the Bianchi identity takes the form 

rfd^H-djP^ = B^G^^n" 
n^d-yP' + \w^djH = B'^'G^^nr (4.16) 

o 

This is no longer symmetric hyperbolic and would not lead to stable constraint preservation even for the Cauchy 
problem. In order to remedy this, the BSSN evolution system modifies (j4.15p by mixing in a set of auxiliary constraints, 
which combine with the constraint system (I4.16P to form a larger symmetric hyperbolic constraint system. The freedom 
in the constraint-mixing parameters and gauge conditions complicates a general treatment. Here the discussion will 
be limited to a particular choice ^38] for which the linearization off Minkowski space leads to a symmetric hyperbolic 
system with a well-posed IB VP for the case of a Dirichlet boundary condition jS'^Ni = on the normal component 
of the shift. Although the nonlinear evolution system is no longer symmetric hyperbolic, the boundary conditions 
for the linearized theory can be formally applied and lead to a symmetric hyperbolic constraint system. Constraint 
preservation then follows for the parameter range {bi < 1, 62 < 1) in the boundary conditions given in equation (97) 
of [11|. The particular choice 61 = 0, leads to the boundary condition fsoj 

H -3P'Ni = G^^n''{n'' ^SN"") ^ Z, (4.17) 

where Z represents contributions from the auxiliary constraints, or, by using the evolution system (|4.15p . 

G^^L^'L'' = Z. (4.18) 

It is a bizarre feature of the 3 + 1 problem that the constraint preserving boundary conditions switch from the outgoing 
Raychaudhuri form (|4.13p to the ingoing Raychaudhuri form (|4.18p in going from the ADM to the BSSN system. The 
Raychaudhuri equation for the outgoing null direction cannot be imposed in the allowed range of (61, 62). Nevertheless, 
(I4.18P can still be used to supply boundary values for the remaining Q^'^k^,y component of extrinsic curvature. 

It appears from the above discussion that the formal application of the Sommerfeld boundary conditions to the 
BSSN system must be restricted to p.l7p . (I3.18P and (|3.20p which supply boundary values for 5 components of the 
extrinsic curvature. Boundary values for the remaining Q^'^fc^^ component must be obtained in accord with constraint 
preservation, e.g. from (j4.18l ) or some variant depending upon the particular formulation. The normal component 
of the shift requires a Dirichlet boundary condition that determines its sign, e.g. (3^Ni — 0. Appropriate boundary 
conditions for the lapse and tangential components of the shift depend upon the specific gauge conditions (see [33| 
for an example). In the spirit of the BSSN formalism, computational experiments would be necessary to determine 
whether p.l7p . (|3.18p and (|3.20p lead to improved performance. 
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V. SUMMARY 

We have shown that a geometricaUy unique spacetime can be locally constructed from Sommerfeld data determined 
by the fields {i, q, ct"'') on a disembodied boundary T, along with the initial data prescribed on Sq- The boundary 
data specify a foliation Bt of T but involve no metric or other geometric properties of the boundary. After the 
embedding oiTuSo as the boundary TuSq of a 4-manifold M, the induced fields {t, q, q°-, a""^) supply the necessary 
boundary data for an isometry class of spacetime metrics which satisfy Einstein's equations. Under diffeomorphisms 
of the disembodied boundary T and Cauchy hypersurface Sq, the mapped data determine diffcomorphic vacuum 
spacetimes. 

The gauge in which a particular metric gat is constructed via a solution of the harmonic IBVP depends upon the 
choice of evolution field t"', which also determines an associated background metric by Lie transport of the initial 
data. Together gab and 17^^ supply the geometric interpretation of the boundary data in terms of the outgoing null 
vector normal to the foliation Bt of the boundary. The field q'^, constructed from q and via (|3.12p . is the 
geodesic curvature of if", relative to its geodesic curvature computed with the background metric. The field a"''' (or 
equivalently a) computed from o'"'' via p.lSp (or (|3.16p ) is the shear of X", relative to its background value. The 
resulting metric is harmonic with respect to the background metric, according to C = (see p.26p ). All possible 
choices of t° are related by diffeomorphism, so that all possible gauges are included. 

The Sommerfeld boundary conditions have direct application to harmonic evolution codes used in simulating binary 
black holes, where they would provide a numerical algorithm based upon a strongly well-posed IBVP. Furthermore, 
for harmonic evolution, homogeneous Sommerfeld data gives rise to asymptotically small back reflection of quadrupole 
waves from an asymptotically large spherical outer boundary of an isolated system. 

The geometric nature of the 6 boundary conditions (|3.17p - (|3.20[) suggests that they might also be applicable 
to codes based upon a 3 + 1 formulation, e.g. the BSSN system, as a way to reduce spurious boundary effects 
generated by the naive Sommerfeld conditions now in practice. There are two caveats. First, strong well-posedness 
of the IBVP does not directly apply to any present 3 + 1 system. Second, the additional Sommerfeld boundary 
conditions (|3.21[) - (|3.23p only guarantee preservation of the Hamiltonian and momentum constraints for harmonic 
(or quasi-harmonic) formulations and would have to be replaced in accord with constraint preservation. In the case 
of the particular BSSN formulation shown in [sll to posses a well-posed IBVP in the linearized approximation, there 
are further complications. Instead of the Sommerfeld condition (j3.19p . a Dirichlet boundary condition is required to 
control the sign of the normal component of the shift. The 5 other Sommerfeld conditions (|3.17p . (|3.18p and (|3.20p 
can be used to supply boundary values for 5 components of the extrinsic curvature. In addition, the boundary values 
of the remaining component of extrinsic curvature can be supplied by a constraint preserving boundary condition, 
but there is no apparent way to do this in a Sommerfeld form. These complications would, at the least, lead to more 
spurious reflection from the boundary than for a harmonic code. The results of this paper can perhaps guide further 
experimentation and investigation towards a modification of the BSSN system that adds to its successful features. 
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